Characterizing and tracking single colloidal particles with video holographic 

microscopy 
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We use digital holographic microscopy and Mie scattering theory to simultaneously characterize 
and track individual colloidal particles. Each holographic snapshot provides enough information to 
measure a colloidal sphere's radius and refractive index to within 1%, and simultaneously to measure 
its three-dimensional position with nanometer in-plane precision and 10 nanometer axial resolution. 



o 

^— > ' 

Or 

q : 
o ■ 

<Z3 ■ 

>v 

^ : 



> ■ 
oo : 

m . 
r- ■ 

(N : 

: 

o ■ 



X 



In addition to their ubiquity in natural and industrial 
processes, colloidal particles have come to be prized as 
building blocks for photonic and optoelectronic devices, 
as probes for biological and macromolecular processes, 
and as model systems for fundamental studies of many- 
body physics. Many of these existing and emerging ap- 
plications would benefit from more effective methods for 
tracking colloidal particles' motions in three dimensions. 
Others require better ways to measure particles' sizes 
and to characterize their optical properties, particularly 
if these measurements can be performed on individual 
particles in situ. 

This Article demonstrates that images obtained with 
in-line holographic microscopy^ can be interpreted with 
Lorenz-Mie theory 3 -! 4 - to obtain exceptionally precise 
measurements of individual colloidal spheres' dimensions 
and optical properties^ while simultaneously tracking 
their three dimensional motions with nanometer-scale 
spatial resolution at video rates^. This method works 
over the entire range of particle sizes and compositions for 
which Mie scattering theory applies, and requires only a 
single calibration of the optical train's magnification. Un- 
like other light scattering techniques for measuring parti- 
cle size^ or refractive index, holographic particle analysis 
can be applied directly to individual particles in hetero- 
geneous samples and also is compatible with scanned 9 
and holographic*^ optical trapping. 

Our holographic analysis instrument is based on a 
standard inverted optical microscope (Nikon TE-2000U), 
with a collimated and attenuated HeNe laser (Uniphase 
5 mW, A = 0.632 /Ltm) replacing the conventional in- 
candescent illuminator and condenser. As indicated 
schematically in Fig.Q] light scattered by a particle prop- 
agates to the microscope's focal plane, where it interferes 
with the undiffracted portion of the beam. The resulting 
interference pattern is magnified 1 by the microscope's ob- 
jective lens (Nikon 100 x NA 1.4 oil immersion Plan-Apo) 
and video eyepiece (1.5 x) onto the sensor of a grey-scale 
video camera (NEC TI-324AII). This system provides a 
total magnification of 135±1 nm/pixel over a 86 x 65 fim 2 
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FIG. 1: Principle of video holographic microscopy. A 

colloidal particle scatters a portion E s (r) of an initially col- 
limated laser beam Eo(r). The scattered beam, here repre- 
sented by 5 calculated iso-amplitude surfaces, interferes with 
the unscattered portion of the incident beam in the focal plane 
of a microscope objective, thereby forming an in-line holo- 
gram, I(p). 



field of view. Images are recorded as uncompressed digi- 
tal video at 30 frames/s using a commercial digital video 
recorder (Pioneer 520HS). 

Analyzing these digitized holograms yields the parti- 
cle's three-dimensional position, r p , its radius, a, and 
its index of refraction, n p . We assume that the incident 
field, Eq(t) = uo{p) exp(ikz) e, is uniformly polarized in 
the e direction and varies slowly enough over the size of 
the particle to be treated as a plane wave propagating 
along the z direction. Its amplitude urj(p) at position 
p = (x, y) in the plane z = z p of the particle is thus the 



same as its amplitude in the focal plane, z = 0. The 
wave propagates along the z direction with wave number 
k = 2im m / \, where A is the light's wavelength in vacuum 
and n m is the refractive index of the medium. For pure 
water at 25°C, n m — 1.3326 at A = 0.632 nm. 

The particle at r p scatters a portion of the incident 
field into a highly structured outgoing wave, E s (r) = 
a exp(-ikzp) u (r p ) f s (r — r p ), where a « 1 accounts 
for variations in the illumination, and where f s (r) is the 
Lorenz-Mie scattering function 3 -^ 1 -, which depends on 
a, n p , n m and A. The scattered field generally covers a 
large enough area at the focal plane that the interference 
pattern, 
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is dominated by long- wavelength variations in |uo(p)| • 
The resulting distortions have been characterized 12 , but 
were not corrected in previous analyses of J(p)£i£i£ i 12 i 13 . 
Fortunately, \uo(p)\ can be measured in an empty field 
of view, and the in-line hologram can be normalized to 
obtain the undistorted image 
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on the plane z = 0. If we further assume that the phase of 
the collimated incident beam varies slowly over the field 
of view, the normalized image is related to the calculated 
Mie scattering pattern, f s (r), in the plane z — by 

B(p) » l+2aM{f s (r - r p ) ■ ee^j+a 2 |/> - r p )\ 2 . 

(4) 

Equation (J4j) can be fit to measured holograms by 
treating the particle's three-dimensional position, its ra- 
dius and its refractive index as free parameters. Previous 
studies fit non-normalized holograms to phenomcnolog- 
ical models ^ 13 ' 14 ' 15 ' 16 ' 17 or Mie scattering theory^ for 
some of these quantities, but never all five. Because 
errors in the adjustable parameters are strongly corre- 
lated, failing to optimize them all simultaneously yields 
inaccurate results. Fitting instead to the full Lorenz- 
Mie theor y 3 ' 4 ' 11 ' 12 ' 19 ' 20 provides more information with 
greater precision. 

Numerical fits to digitized and normalized holographic 
images were performed with the Levenberg-Marquardt 
nonlinear least-squares minimization algorithm 21 i 22 i 23 
using the camera's measured signal-to-noise ratio to esti- 
mate single-pixel errors. The \ 2 deviates for all of the fits 
we report are of order unity, so that the calculated un- 
certainties in the fit parameters accurately reflect their 
precisio n 21 ' 23 ' 24 . These estimates incorporate the esti- 
mated covariance of the adjustable parameters, so that 
they also may be interpreted as the resolution of each 
parameter 
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FIG. 2: Fitting to normalized holograms, (a) Nor- 
malized hologram B{p), numerical fit to Eq. Q, and az- 
imuthally averaged radial profile B(p) for a 1.43 fim diameter 
polystyrene sphere in water at z p = 22.7 /im. All scale bars 
indicate 10 /xm. Curves in the radial profile are obtained from 
experimental data, discrete points were obtained from the fit. 
(b) Data for a 1.45 jj,m diameter Ti02 sphere dispersed in 
immersion oil (n m = 1.515) at z v = 7.0 /im (c) Data for a 
4.5 fj,m diameter SiC>2 sphere in water at z p = 38.8 fim. 



Because the laser's wavelength and the medium's re- 
fractive index are both known, the only instrumen- 
tal calibration is the overall magnification. This con- 
trasts with other three-dimensional particle tracking 
techniques 1 -' 2 - ! 13 ' 2 ^ 26 ' 27 , which require independent cali- 
brations for each type of particle, particularly to track 
particles in depth. 

The image in Fig. [Ha) shows the normalized hologram, 
B(p), for a polystyrene sulfate sphere dispersed in wa- 
ter at height z p = 22.7 ^m above the focal plane. This 
sphere was obtained from a commercial sample with a 
nominal diameter of 2a = 1.48 ± 0.03 /im (Bangs Labs, 
Lot PS04N/6064). The camera's electronic shutter was 
set for an exposure time of 0.25 msec to minimize blur- 
ring due to Brownian motion^. After normalizing the 
raw 8-bit digitized images, each pixel contains roughly 5 
significant bits of information. The numerical fit to B(p) 
faithfully reproduces not just the position of the inter- 
ference fringes, but also their magnitudes. The quality 
of the fit may be judged from the azimuthal average; 
the solid curve is an angular average about the center of 
B(p), the dashed curves indicate the standard deviations 
of the average, and the discrete points are obtained from 
the fit. 

The fit value for the radius, a — 0.73 ± 0.01 p,m, falls 
in the sample's specified range, which reflected a lower 
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FIG. 3: Holographic tracking of a sedimenting colloidal silica sphere, (a) and (c) DHM images of the sphere at 
the beginning and end of its trajectory, respectively. The scale bar indicates 5 /im. (b) and (d) Fits to Eq. (e) Three- 
dimensional trajectory with starting point (circle) and end point (square) labeled, (f) z(t), showing thermal fluctuations about 
uniform sedimentation. Inset: The fit refractive index is independent of position, (g) Mean-square positional fluctuations 
display Einstein-Smoluchowsky scaling in x, y and z. 



bound of 0.69 ± 0.07 /um obtained with a Beckman Z2 
Coulter Counter and an upper bound of 0.76 ± 0.08 
obtained by analytical centrifugation. Agreement be- 
tween the quoted and measured particle size suggests 
that the present measurement's accuracy is comparable 
to its precision. In that case, both precision and accuracy 
surpass results previously obtained 6 through analysis of 
I{p). The trajectory-averaged value for the refractive 
index, n p = 1.55 ± 0.03, also is consistent with the prop- 
erties of polystyrene colloid inferred from light scattering 
measurements on bulk dispersions^ 9 -. 

Comparable precision in measuring a single particle's 
refractive index has been achieved by analyzing a col- 
loidal particle's dynamics in an optical trap 30 . This 
method only can be applied to particles with compara- 
tively small refractive indexes, however, because particles 



with relative refractive indexes greater than n p ss 1.3 n m 
are difficult to trap. Holographic characterization, by 
contrast, requires only a single holographic snapshot 
rather than an extensive time series, does not require 
optical trapping, and so does not require separate cali- 
bration of the trap, and is effective over a wider range of 
particle sizes and refractive indexes. 

The corresponding data in Fig. [2fb) were obtained for 
a 1.45 /im diameter Ti02 sphere at z p = 7 /jm above 
the focal plane. This sample was synthesized from ti- 
tanium tetraethoxide and was heat-treated to increase 
its density 3 ^. Strong forward scattering by such high- 
index particles gives rise to imaging artifacts unless the 
medium is index matched to the cover slip. Dispersing 
the particle in immersion oil [n m = 1.515) eliminates 
these artifacts, but introduces spherical aberration for 
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the lens we used, which must be corrected", to obtain 
reliable results. The fit diameter of 1.45 ± 0.03 /*m and 
refractive index of 2.01 ± 0.05 are consistent with results 
obtained by electron microscopy and bulk light scatter- 
ing, respectively. This result is noteworthy because no 
other single-particle characterization method works for 
such high refractive indexes. 

The data in Fig.^c) show results for a nominally 5 /im 
silica sphere (Bangs Labs, Lot SS05N/4364) dispersed in 
water at z p — 38.8 /zm above the focal plane. The fit 
refractive index, n p = 1.434 ± 0.001, is appropriate for 
porous silica and the diameter, a = 4.51 ±0.01 pm agrees 
with the 4.82 ± 0.59 /im value obtained for this sample 
with a Beckman Z2 Coulter Counter. 

The same fits resolve the particle's position with a pre- 
cision of 1 nm in-plane and 10 nm along the optical axis. 
This substantially improves upon the typical 10 nm in- 
plane accuracy obtained with standard particle tracking 
techniques with the same microscope and camera^ 5 -. The 
difference can be ascribed to the larger number of pixels 
subtended by a holographic image, and to the images' 
strong intensity gradients, which constrain the fits. The 
estimated 10 nm axial resolution surpasses results ob- 
tained with morphometric axial particle trackin g 1 ' 2 ! 25 by 
a factor of ten. 

Nanometer-scale tracking resolution can be obtained 
under conventional illumination, but requires detailed 
calibrations for each particle^. Still better in- 
plane spatial resolution can be obtained at much 
higher bandwidths through back-focal-plane interfero- 
metric methods^, but also require accurate calibrations 
with piezo translators. Total internal reflection mi- 
croscopy (TIRM) similarly offers sub-nanometer axial 
resolutio n 35 ' 36 , but performs no better than conventional 
imaging methods for in-plane tracking. 

An additional benefit of holographic imaging over 
other particle-tracking techniques is its very large depth 
of focus. Our system provides useful data over a range 
of more than 100 /im, which contrasts with the ±3 /nm 
useful depth of focus using conventional illumination 33 
and the 100 nm range of TIRM 3 ^. 

Holographic video microscopy lends itself to three- 
dimensional particle tracking, as the data in Fig. [3] 
demonstrate for a colloidal silica sphere (Bangs Labs, Lot 
SS04N/5252) dispersed in water. This particle was lifted 
30 fim above the focal plane with an optical tweezer, and 
then released and allowed to sediment. The images in 
Fig. E£a) and (c) show the particle near the beginning of 
its trajectory and near the end. Fits to Eq. ([4]) are shown 
in Figs. [Sib) and (d). 

The particle's measured trajectory in 1/30 s intervals 
during 15 s of its descent is plotted in Fig.[3]Je). Its verti- 
cal position z(t), Fig.[3ff), displays fluctuations about a 
uniform sedimentation speed, v — 1.021 ± 0.005 pm/s. 
This provides an estimate for the particle's density 
through p p = p m + 9riv/ (2a 2 g), where p m = 0.997 g/cm 3 
is the density of water and rj = 0.0105 P is its viscosity 
at T = 21°C, and where g = 9.8 m/s 2 is the acceleration 



due to gravity. The fit value for the particle's radius, at 
a = 0.729 ± 0.012 /im, remained constant as the parti- 
cle settled. This value is consistent with the manufac- 
turer's specified radius of 0.76 ± 0.04 /im, measured with 
a Beckman Z2 Coulter Counter. Accordingly, we obtain 
p p = 1.92 ± 0.02 g/cm 3 , which is a few percent smaller 
than the manufacturer's rating for the sample. However, 
the fit value for the refractive index, n p = 1.430 ± 0.007, 
also is 1.5% below the rated value, suggesting that the 
particle is indeed less dense than specified. 

The mean-square displacements, Ar 2 (r) = 

(( r j(t + T ) — r jW) 2 )> °f the components of the 
particle's position provide additional consistency 
checks. As the data in Fig. GUg) show, fluctuations 
in the trajectory's individual Cartesian components 
agree with each other, and all three display linear 
Einstein-Smoluchowsky scaling, Ar 2 (r) = 2Dt, with 
a diffusion coefficient D = 0.33 ± 0.03 /im 2 /s. This is 
consistent with the anticipated Stokes-Einstein value, 
D = k B T/{6nr]a) = 0.30 ± 0.02 pm 2 /s, where k B is 
Boltzmann's constant. Using the methods of Refi 2 ^, we 
then interpret the offsets obtained from linear fits to 
Ar 2 (i) to be consistent with no worse than 1 nm accu- 
racy for in-plane positions and 10 nm for axial positions 
throughout the trajectory. The optical characterization 
of the particle's properties thus is consistent with the 
particle's measured dynamics. 

We have successfully applied holographic characteriza- 
tion to colloidal spheres as small as 100 nm in diameter 
and as large as 10 /jm. Unlike model-based analytical 
methods, fitting to the exact Lorenz-Mie scattering the- 
ory is robust and reliable over a far wider range of particle 
sizes, provided that care is taken to maintain numerical 
stability in calculating / s (r)£ i 19 ' 20 . Such numerical im- 
plementations have been reported for particles as small as 
a few nanometers and as large as a few millimeters, with 
relative refractive index ratios from less than m = 0.1 
to over 10, and with large imaginary refractive indexes. 
In all cases, the instrumental magnification and field of 
view must be selected to fit the sample. 

The principal limitations of the six-parameter model in 
Eq. (TJ| are the assumptions that the scatterer is homoge- 
neous and isotropic, and that its interface is sharp. These 
assumptions can be relaxed at the cost of increased com- 
plexity and reduced numerical robustness. For example, 
analytical results are available for core-shell particles 3 , 
and for particles with more complex shapes&ii, such as 
ellipsoids, spherical clusters and cylindrical nanowires. 
All such elaborations involve additional adjustable pa- 
rameters and thus are likely to pose computational chal- 
lenges. 

We have demonstrated that a single snapshot from an 
in-line holographic microscope can be used to measure 
a colloidal sphere's position and size with nanometer- 
scale resolution, and its refractive index with precision 
typically surpassing 1 percent. 

A video stream of such images therefore constitutes a 
powerful six-dimensional microscopy for soft-matter and 
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biological systems. Holographic particle tracking is ideal 
for three-dimensional microrheology, for measuring col- 
loidal interactions and as force probes for biophysics. 
The methods we have described can be applied to track- 
ing large numbers of particles in the field of view simul- 
taneously for highly parallel measurements. Real-time 
single-particle characterization and tracking of large par- 
ticle ensembles will be invaluable in such applications 
as holographic assembly of photonic device a 37 i 38 . Ap- 
plied to more highly structured samples such as biolog- 
ical cells and colloidal heterostructures, they could be 



used as a basis for cytometric analysis or combinatorial 
synthesis 3 ^. 
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